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In this work, we present an analytical theory of strongly correlated Wigner crystals (WCs) in 
the lowest Landau level (LLL) by constructing an approximate, but accurate effective two-body 
interaction for composite fermions (CFs) participating in the WCs. This requires integrating out 
the degrees of freedom of all surrounding CFs, which we accomplish analytically by approximating 
their wave functions by delta functions. This method produces energies of various strongly correlated 
WCs that are in excellent agreement with those obtained from the Monte Carlo simulation of the 
full CF crystal wave functions. We compute the compressibility of the strongly correlated WCs 
in the LLL and predict discontinuous changes at the phase boundaries separating different crystal 
phases. 

PACS numbers: 73.43.-f, 73.50.-h, 71.10.Pm 


Wigner predicted more than eighty years ago [1] that, 
when the interaction energy is dominant over the ki¬ 
netic energy, electrons form a crystal, which is called 
WC after its originator. One possible way of suppress¬ 
ing the kinetic energy relative to the interaction energy 
is via the application of a strong magnetic field to two- 
dimensional electron systems [2], which generates a fas¬ 
cinating series of various emergent quantum phases. The 
most celebrated examples are the fractional quantum 
Hall states [3], where new emergent quasiparticles called 
CFs form the quantum Hall liquid states [4, 5]. The 
quantum Hall liquid states are more effective in minimiz¬ 
ing the interaction energy than WCs for a range of filling 
factors that is not too low. Nevertheless, WCs are ex¬ 
pected to occur at sufficiently low filling factors. Indeed, 
insulating states observed at filling factor z/ < 1/5 are 
interpreted as pinned WCs [6-20]. More recently, indi¬ 
cations of the existence of a WC in the LLL have been 
seen through commensurability magneto-resistance oscil¬ 
lations in bilayer Hall systems composed of a CF sea in 
one layer and a WC in the other [21]. 

Numerous theoretical studies have investigated the na¬ 
ture of WCs in the LLL [22-37]. Initially, Maki and 
Zotos [22] considered an uncorrelated Hartree-Fock WC 
of electrons, which was improved upon by Lam and 
Girvin [23] by incorporating correlations. In view of the 
success of the CF theory, Yi and Fertig [27] proposed 
a strongly correlated WC composed of CFs, which was 
subsequently shown by Chang et al [33] to provide an 
accurate description at low filling factors (z/ < 1/5). De¬ 
spite these extensive theoretical works, the calculation of 
a precise phase diagram of quantum Hall liquids versus 
CF crystals (CFCs) remained stalled for many years due 
to difficulties in obtaining the energy of CFCs in the ther¬ 
modynamic limit accurately. This issue was resolved in a 
recent work [38] inspired by the Thomson problem [39]. 
Here, the CFG wave functions are constructed in the 
spherical geometry by placing the WC wave packet cen¬ 


ters at the locations that minimize the Coulomb energy 
of N charged point particles on the surface of a sphere. 
Locally, these minimum energy positions resemble the 
hexagonal lattice, which is the minimum energy symme¬ 
try for a classical 2D electron crystal [40]. This allows a 
precise investigation of the CFG wave functions up to a 
fairly large system size (A^ ^ 100) [38]. 

The Monte Carlo (MC) simulation of the CFC wave 
functions is computationally quite expensive and rather 
difficult to implement. Furthermore, it turns out that 
even though the energy obtained from this method en¬ 
ables a determination of the phase diagram, it is not suffi¬ 
ciently accurate to allow an evaluation of quantities such 
as compressibility, which is related to the second deriva¬ 
tive of the energy. We develop in this work an analytical 
theory of the CFCs by constructing an accurate effective 
two-body interaction, which is based on the two-body 
wave function of CFs participating in the CFCs. This 
requires integrating out the degrees of freedom of all 
surrounding CFs, which we accomplish analytically by 
approximating their wave functions as delta functions. 
We call this approach the “renormalized two-body for¬ 
malism,” to be contrasted with the “isolated two-body 
formalism” where the effects of all surrounding CFs are 
neglected. The CFC energies obtained from the renor¬ 
malized two-body formalism are in excellent agreement 
with those obtained from the MC simulation of the full 
CFC wave functions. With these analytical results, we 
obtain the compressibility and predict that its measure¬ 
ments, such as those carried out in GaAs heterostructures 
or in graphene [41-46], can detect the phase diagram of 
the CFCs. 

We begin by constructing the wave function for the 
CFC carrying 2p vortices, or in short ^^CFC, as follows: 

= n(% - (1) 

j<k 

where Zj = Xj + ipj is the coordinates of the j-th elec¬ 
tron. z/ and z/* denote the filling factors of electron and 
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CFCs, respectively. The function = Det [(^r. (r^)] 
is the MZ wave function for the uncorrelated WC, com¬ 
prising the LLL coherent-state wave function ^R(r) = 
exp [—(r — R)^/4 — i[xY — ^X)/2] centered at R = 
(X>) [22]. 

Without the Jastrow factor correlation, the energy of 
the uncorrelated Maki-Zotos (MZ) WC, i.e., ^CFC, can 
be computed from the effective two-body interaction [22]: 


e^/e/e 




|ri-r2|-i 



= ^sech(i?y8)/o(i?y8), 


( 2 ) 


where Rij = |Rij| = jR^ — Rj| is the distance be¬ 
tween two crystal lattice centers. The uncorrelated 
two-body wave function is given by ^^^(i*i,r 2 ) = 
Co-4.[0r. (ri)0R^. (r 2 )] with Co being the normalization 
constant and A being the antisymmetrization operator. 
Jo is the modified Bessel function of the first kind. The 
energy per particle of the uncorrelated MZ WC, is 

computed by performing the Madelung-type lattice sum¬ 
mation of the MZ effective two-body interaction energy 
between all pairs of electrons in the hexagonal lattice: 



^MZ 

e^/e/e 


1 1\ 
2^^V eyelB Rij) 


ay/v, (3) 


where the terms 2 ^ + cty/v with a = 0.782133 

are subtracted to take into account the neutralizing effect 
of the uniform positive-charge background [40]. 

For strongly correlated WCs, we need to take care 
of the Jastrow factor. As a first try, we begin by 
focusing on two isolated CFs and ignoring all other 
surrounding CFs, in which situation the two-body 
CF wave function is obtained as r 2 ) = 

C 2 p{zi - 2 : 2 )^^^ [0R,(ri)0R^.(r2)]. The normaliza¬ 
tion constant is C 2 p = [iFi(2p -h l;l;Rf^/4) — 

iFi(2p+ 1; 1; -R%/4)]-^/^^2V{2p 
where iFi(a; b; z) = i® Rummer’s hy¬ 

pergeometric function with = a(a-l-l) • • • (a-|-n —1). 
We refer to this approach the isolated two-body formal¬ 
ism. Figure 1 (a) shows a schematic diagram for the iso¬ 
lated two-body formalism, which is accompanied by the 
probability density of the two-body CF wave function 
defined by Pisoi(r) = in Fig. 1 (c). 

It is important to note that the actual distance be¬ 
tween composite fermions dij is not exactly equal to the 
nominal distance Rij^ because the Jastrow factor incor¬ 
porates an additional repulsion into the two-body wave 
function, pushing the wave packets slightly farther apart. 
To take this effect into account, we define dij to be the 
median distance between two maxima in the two-body 
probability density. Specifically, dij is related with Rij 
so that = 0 at r = ri - r 2 = (dy,0) 


FIG. 1: (Color online) Schematic diagram showing (a) the iso¬ 
lated and (b) the renormalized two-body formalism, accompa¬ 
nied by the probability density of the two-body CF wave func¬ 
tion in each formalism, (c) pisoi(r) = / 

and (d) prenorm(r) = / d'^v', T ) f, rcspectivcly. 
Here, we set 2p — 2 and the nominal distance between crystal 
centers to be 6/b. Note that probability densities are plotted 
in the natural log scale. 


with Hij = (Rij,0), which, after some algebra, becomes 
d‘fj — Rij dij coth {Rij dij/A) — 8p = 0. Note that dij is 
slightly larger than Rij with their difference growing as 
2p increases. Given this information, the effective two- 
body interaction between CFs in the isolated two-body 
formalism is computed as follows: 


^isol ) 

e^/ c/b 







= B. 


2p- 


iFi(2p + 1/2; 1; R^/A) - ifi(2p -h 1/2; 1; -R^/A) 


L2p{-RlJ4)e^y^ - L2p(i?? /4)e-««^ 


( 4 ) 


where B 2 p = r(2p H- l/2)/[2r(2p+ 1)] and Ln{x) is the 
Laguerre polynomial. For the hexagonal lattice, dij is 
set equal to the distance between various crystal lattice 
centers via dij = + pb‘^ with a = {Yk j 

and b = y/3a. 

In the isolated two-body formalism, the energy per par¬ 
ticle of ^^CFC, evaluated similarly to Eq. (3) 

by replacing with Figure 2 (a) shows 
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FIG. 2: (Color online) Energy per particle of various CFG 
states, E ^ in reference to that of the MZ WC, as a 

function of filling factor, z/, obtained in (a) the isolated and (b) 
the renormalized two-body formalism. Symbols and dashed 
lines are obtained from the MG simulation of the full CFG 
wave functions. Note that each ^^CFC is the lowest energy 
state (at least, among various CFG states) within the range of 
l/(2p+3) <iy< l/(2p+l) for 2p > 2. The MZ WC, or °CFC 
is not energetically favorable at any filling factor ranges. 


as a function of filling factor, which is com¬ 
pared with the MC simulation results obtained from the 
full CF wave function in the spherical geometry. As one 
can see, — E^^ shows a reasonably good agree¬ 

ment with the MC simulation results, especially at low 
filling factors. There are, however, some sizable quanti¬ 
tative discrepancies in general. 

To improve upon the isolated two-body formalism, it 
is necessary to include Jastrow-factor correlation effects 
arising from all surrounding CFs in some fashion. To 
this end, we test an approximation called the surround¬ 
ing delta-function approximation^ where the wave func¬ 
tions of all surrounding CFs are approximated as delta 
functions. See Fig. 1 (b) for a schematic diagram. This 
approximation should be exact in the limit of large sepa¬ 
ration between crystalline CFs. Within the surrounding 
delta-function approximation, one can integrate out the 
degrees of freedom of all surrounding CFs and then derive 
the analytical two-body CF wave function. For conve¬ 
nience, we call this approach the renormalized two-body 
formalism. 

In the renormalized two-body formalism, the two- 
body CF wave function is written as '0reS)^m(^i5 ^ 2 ) cx: 

r 2 ) ni^i,jiz 2 -ZirP , where 

Zk denotes the coordinates of the ^-th crystal lattice 


center. After dividing a constant factor Ylf.^ij{Zi — 
Zk)‘^P Yli^. j{Zj — Zi)‘^P, The above equation can be 
rewritten as 


V’reforS(ri,r2) = C2p(^l “ ^2)^^ 


X A 


r-f(ri,r2)(/R,(ri)</R (r2) 


( 5 ) 


where C 2 p is the normalization constant, r^J(ri,r 2 ) = 
{Zi - Z^rp/[{z, - Z,) 2 j >(^2 - Zi)^% and = 

{E-zI^p - ‘/'R.W [.^R.(r)]^, which is the 
renormalized version of the coherent-state wave function 
centered at R^. 

It is important to note that all the complicated many- 
body correlations are embedded in the renormalization 
factor Er. (r). The success of this approach stems from 
the fact that we are able to obtain the analytical form of 

jr e^{l{z-z,)\i^Je,{f{z-z,)\ip 

^ Js{liz-z,)\r-)es{f{z-z,)\ip 

^3(o|i!)03(o|if) ’ 


where On{z\T) is the Jacobi theta function and b = 
VSa with a being the lattice constant. See Supple¬ 
mentary Material for details. Figure 1 (d) shows the 
renormalized probability density defined by Prenorm(i*) = 
J . As one can see, the renormalized 

probability density exhibits the hexagonal symmetry of 
the WC which was absent in the isolated two-body for¬ 
malism. 

As before, the energy per particle in the renormal¬ 
ized two-body formalism, can be computed as 

the Madelung-type lattice summation of the renormal¬ 
ized effective two-body interaction, Kenorm/(^^/^^ b) = 
i^TpormW^i - ^otc that, in general, 

Kenorm depends on the vector not just on the 

distance Rij. Conveniently, however, f^enorm can be 
well approximated as a function of only Rij if Rij > 
[47r(2p-b 1 )/a/ 3]^/^. Also, due to the additional repul¬ 
sion from all surrounding CFs, the actual distance dij 
becomes quite close to the nominal distance Rij in most 
situations so that Rij can be simply regarded as dij. 

Figure 2 (b) shows — E^^ as a function of 

filling factor. As one can see, the results from the renor¬ 
malized two-body formalism are in excellent agreement 
with those from the MC simulation of the full CFC wave 
functions. It is interesting to observe that the renor¬ 
malized two-body formalism can even capture the initial 
upturn of the energy near n = l/{2p +1) for each cor¬ 
responding ^^CFC. The most significant discrepancy is 
that the MC results exhibit sharp drops immediately fol¬ 
lowing such upturns. It is important to note, however, 
that, regardless of being isolated or renormalized, the 
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FIG. 3: (Color online) Shear modulus, (7^, of various CFG 
states as a function of filling factor, z/, obtained in (a) the 
isolated and (b) the renormalized two-body formalism. For 
convenience, is expressed as its relative ratio with respect 
to that for the classical WC, /els) — 0.0978v^. 

Note that of each ^^CFC is valid only within the range of 
l/(2p + 3) < < l/(2p + 1), where its curve is plotted in a 

solid line, denoting that the ^^CFC is the lowest energy state 
here. 


two-body formalism for each ^^CFC is supposed to lose 
its validity near v = l/(2p+ 1) since, here, wave packets 
are highly overlapping and thus higher-body corrections 
become important. Given the simplification in the two- 
body formalism, we consider the agreement to be excel¬ 
lent. 

Bolstered by the quantitative accuracy of the renor¬ 
malized two-body formalism, we now compute the shear 
modulus of CFG states, as a function of filling fac¬ 
tor. To this end, we utilize the following relation [38]: 

is important to note that this 
relation is derived under the assumption that only the 
two-body interaction is relevant, and is therefore consis¬ 
tent with our two-body formalism. Figure 3 shows of 
various CFG states as a function of filling factor obtained 
by using and While obtained from 

the isolated two-body formalism is not to be trusted, it 
shows the relative importance of correlations from all sur¬ 
rounding GFs, which enhance significantly. In partic¬ 
ular, in the renormalized two-body formalism shows 
a series of huge enhancements followed by discontinuous 
drops near u = l/{2p +1), which can be used as a dis¬ 
tinctive signature for a phase transition between different 
GFG states. 

Another important observable is the compressibility, 
whose inverse can be computed as follows [47]: = 


FIG. 4: (Golor online) Inverse of the compressibility, z^~^, of 
various GFG states as a function of filling factor, z/, obtained 
in (a) the isolated and (b) the renormalized two-body formal¬ 
ism. Similar to G^, of each ^^GFG is valid only within 
the range of l/(2p + 3) < z/ < l/(2p + 1), where its curve is 
plotted in a solid line. 

where it is used that the electron 

density is related with the filling factor via n = v 
Figure 4 shows of various GFG states as a function of 
filling factor. At first sight, it may seem strange that the 
compressibility becomes negative in some regimes. This 
does not, however, mean an instability here since, by con¬ 
struction, we do not allow the positive background charge 
to relax. What we obtain above is the electronic part of 
the compressibility called the proper compressibility [48], 
which can be negative. In fact, the proper compressibil¬ 
ity is directly measured in capacitive experiments [41-43] 
or by scanning single-electron transistor [44-46]. Gom- 
pressibility has served as a powerful tool for detecting 
phase transitions between different fractional quantum 
Hall states as well as between differently spin polarized 
states at a given fraction [44-46]. Our calculations pre¬ 
dict discontinuous changes in compressibility at the phase 
boundaries separating the different GFG phases, which 
can allow a determination of the phase diagram. Ob¬ 
servation of such transitions inside the crystal phase will 
serve as direct evidence for the correlated GF character 
of WGs in the LLL, corroborating existing experimental 
indications of the GFG states [49, 50]. 

It is noteworthy that this is the first example where an 
accurate analytical treatment has become possible for a 
strongly correlated state in the LLL. It would be worth 
investigating if our method can be extended to the liquid 
states of GFs. 
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Supplementary Material: Theta-function Representation of the Renormalization Factor 


The goal of this Supplementary Material is to evaluate the renormalization 
follows: 


factor (r), which is defined as 




TT {Z - Zkf 

11 rz. - 7.,X‘ 


( 7 ) 


where the product index, k, goes over all crystal lattice points on the triangular lattice except the z-th. Now, for a 
reason that becomes clear below, it is convenient to decompose the triangular lattice into two overlapping rectangular 
lattices, say, the rectangular lattice 1 (RLi) and 2 (RL 2 ). See Fig. 5 for illustration. 

By using this decomposition, Eq. (7) can be rewritten as follows: 

-^R.(r)=-F^^(r)J-Sfnr), (8) 


where (r) and (r) are the renormalization factors containing contributions from RLi and RL 2 , respectively: 

{Z - Zk? 


T 


•RL, 

R 7 


(r)= n 

/e/i,eRLo 


{Zi-ZkY^ 


( 9 ) 


where ci = 1 or 2. Specifically, the coordinates of crystal lattice points can be specified via = ma + inh for 

RLi and = (m — l/2)a + i(n — 1/2)6 for RL 2 , where m and n are integers, and a and b (= VSa) are the lattice 

constants of the rectangular lattices along the x and y direction, respectively. In this notation, Eq. (9) can be written 
as follows: 

00 

•^Rb(r)= n 

m,n= — oo 


1 - 




^RLa 
^ {m,n) 


( 10 ) 


which can be further simplified as follows by redefining (m, n) such that the coordinates of all crystal lattice points 
in the product are measured with reference to R^: 






CX) 

n 

m ^n= — 00 
(m,rz)#(0,0) 


y 

ma + inh J 


00 


n 

m,n= —00 


(i -£—- 

\ (m — l/2)a + i{n — 1/2)6 


2 


(11) 

( 12 ) 



Triangular Lattice 




Rectangular Lattice 1 Rectangular Lattice 2 


FIG. 5: (Color online) Decomposition of the triangular lattice into two overlapping rectangular lattices, say, the rectangular 
lattice 1 (RLi) and 2 (RL 2 ). Note that 6 = VSa. 
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The concrete evaluation of the renormalization factor begins with the following, curious representation of the Jacobi 
theta function in terms of the double infinite product [51]: 


N 


M 


Oi{z\r) = z0[{0\t) lim IT lim IT (l + 


N^oo M^oo 

n=-N m=-M 

(n ,m)^(0,0) 


{m + nr)7r 


N 


M 


03 { z \ t ) = 6»3(0|t) lim JJ lim 1 + 

N^oo M^oo '*--*■ \ 

n=l — N m=l — M 


[m - i + (n - i) r] TT 


(13) 


(14) 


where the order of the limits cannot be switched since the double products above are not absolutely convergent. Note 
that the Jacobi theta functions, Oi{z\r) and 6^3(2:|r), are defined as follows: 


^) = 2 ^(—sin (( 2 n + 1 ) 2 :), 

n=0 

00 

03 {z\t) = 63 {z,q) = 1 + 2 "^ cos { 2 nz), 


(15) 

(16) 


where q = 

As one can see, Eqs. (11) and (12) are very similar to Eqs. (13) and (14), respectively, which suggests the following, 
possible equalities: 


N 


M 


= ihn n j“L n b 


N^oo M^oo 

n=—N m=-M 


(m + mf)7r 




(17) 


and 


N 


M 




lim 

N^oo M^oo 

n=—N m=—M 


n 


':{z-Zi) 


[TO-i + (n-i)i^]7r 



leyi(z-z,)\i^y 

) ■ 

[ ^3 (0|y) J 


(18) 


where dummy product indices, m and n, are reversed in sign and shifted by constant integers when necessary. Despite 
the obvious similarity, however, there is a serious problem, which makes the above equalities ill-defined, or at least 
incorrect in the context that we intend. The problem stems from the previously-mentioned fact that the double 
products in the representation of the Jacobi theta function are not absolutely convergent and thus the limits should 
be taken with care. Technically, Eqs. (17) and (18) can be regarded as being obtained by first multiplying the entire 
Jastrow factors contributed by the composite fermions located along the x direction at a given y coordinate and then 
repeating the same procedure for different y coordinates to cover the whole two-dimensional lattice. Eor convenience, 
let us call this the x-priority limit procedure. 

Another way of covering the whole two-dimensional lattice is to multiply the Jastrow factors along the y direction 
at a given x coordinate first and then repeating the same procedure for different x coordinates. This alternative 
covering can be implemented via a slight rearrangement of the product, generating the following, possible equalities: 


N 


M 


n ii-L n 11 

n=-N 


'.= -M 
n)7^(0,0) 


(^ - Zi) 


(n -h mi|) TT 


y 

- eyfiz - Zi)\iy ' 

1 - 

e[{o\if)f{z-z,)_ 


(19) 


and 


N 


M 


•^f-^rior,R.(r)= lim n Jim n b 


(^ - Zi) 


N^oo M^oo 

n=—N m=—M 


[n-i + (m-i)if]7r 


^ 0,{f(z-Zi)\ify ^ 

^3 (oNf) 


( 20 ) 


where, as before, dummy product indices, m and n, are reversed in sign and shifted by constant integers when 
necessary. In contrast to the x-priority limit procedure, we call this the y-priority limit procedure. It is intriguing to 
check numerically that Eqs. (19) and (20) are actually different from Eqs. (17) and (18), respectively, which presents 
us a dilemma regarding what form to choose. 
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FIG. 6: (Color online) Numerical tests for the validity of the geometric mean conjecture for the renormalization factor, 
^?D.R,(r). Here, is set to be at the origin, i. e., Ri = 0, without loss of generality. We plot in the upper 

panels and |*F^^,^Q(r)| in the bottom panels. As one can see, the geometric mean conjecture for the renormalization factor 
yields essentially exact results in a wide range of r = (x, y) for various values of the lattice constants, a and h. 


In fact, the most natural physical limit is neither the x nor ^-priority limit, but rather to evaluate the Jastrow-factor 
product for a given finite two-dimensional (2D) cluster with the linear size S and increase S while keeping the aspect 
ratio of the cluster. That is to say. 


Ns 


Ms 


^2d‘',r.m=, 1m n n (1- 


n= — Ns rn=-Mg 

Ns Ms / 


=,im n n 


1 - 


z- Zj 
ma + inb , 

^ - Zi 


n=—Ns m= — Ms 


(m - ^) a + i (n - I) 6 


( 21 ) 

( 22 ) 


where Ms = S and Ns = round(|S'). Let us call this the proper 2D limit procedure. 

Now, a big question is how to actually evaluate the renormalization factor in this proper 2D limit. To this end, we 
make a simple conjecture, which we call the geometric mean conjecture. Specifically, considering that the x-priority 
limit [in Eqs. (17) and (18)] and the ^-priority limit [in (19) and (20)] are the two opposite limits with reversed roles 
between the x and y direction, we make a conjecture that the proper 2D limit may be obtained by simply taking the 
geometric mean between the x and ^-priority limits: 


which results in 




I^RLq 


a: —prior,Fli V / y—prior,v 


■^RLi 
*^20,R 




m){z-z,Y 


and 


(23) 


(24) 


^3(f(^ 


^3(01*1)03 (0|*f) 




( 25 ) 
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where we take the liberty of ignoring the phase-factor issue, which may occur when taking the square root. We test 
the validity of the geometric mean conjecture by comparing Eqs. (24) and (25) with the numerical results evaluated 
according to Eqs. (21) and (22), respectively, for sufficiently large S. Eigure 6 shows the comparison between the 
geometric mean conjecture and the numerical results in a wide range of r for various values of the lattice constants, 
a and b. As one can see, the agreement is essentially exact within numerical accuracy, establishing the validity of the 
geometric mean conjecture. 

In conclusion, we obtain the final form of the renormalization factor as follows: 


7-R.(r) = (r) = 


RL2 


- zm) (f - zm) e, - z,)\i^) e, (f ) 
(ON!) (ON!) 0 , (o|i^) 03 (o|if) 


(26) 
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